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Resumen 


La conservación del suelo y el agua es una prioridad en las cuencas 
hidrográficas de las regiones áridas y semiáridas para la adecuada 


planificación y gestión integrada de los recursos hídricos. El objetivo de 
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este trabajo fue priorizar 91 cuencas hidrográficas en 14 regiones del Perú 
con diferentes condiciones geográficas, hidrológicas y geológicas, 
mediante la integración de sistemas de información geográfica (SIG), 
análisis de componentes principales (PCA) y el análisis de suma 
ponderada (WSA). Se identificaron regiones homogéneas con el análisis 
de conglomerados jerárquico en el método de R y Ward. El resultado 
mostró la existencia de un 19.49 % del área total en categoría de alta y 
muy alta prioridad en dos regiones. Los resultados del análisis de 
conglomerados mostraron que el 35 % de las cuencas hidrográficas tiene 
zonas homogéneas dentro del factor geométrico y de forma, mientras que 
el 65 % dentro del factor de drenaje y relieve. En general, los métodos 
SIG, PCA y SWA son herramientas eficientes, que permiten a las 
autoridades encargadas la toma de decisiones para una mejor 
planificación y conservación de los recursos naturales en las cuencas 


hidrográficas. 


Palabras clave: erosión, morfometría, drenaje, análisis de 


conglomerados, análisis factorial. 


Abstract 


Soil and water conservation is a priority in the watersheds of arid and 
semi-arid regions for the proper planning and integrated management of 
water resources. The objective of this work is prioritizing 91 watersheds 
in 14 regions of Peru with different geographical, hydrological and 
geological conditions, through the integration of Geographic Information 
Systems (GIS), Principal Component Analysis (PCA) and the Weighted 


Sum Approach (WSA). In addition, homogeneous regions were identified 
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with hierarchical Cluster Analysis in R and Ward's method. The result 
showed the existence of 19.49 % of total area in high and very high 
priority category in two regions. The results of the Cluster Analysis 
showed that 35 % of the watersheds have homogeneous zones within the 
geometric and shape factor, while 65 % within the drainage and relief 
factor. In general, GIS, PCA and SWA methods are an efficient tool, which 
allows decision-making authorities for better planning and conservation 


of natural resources in hydrographic watersheds. 


Keywords: Erosion, morphometric, drainage, cluster analysis, factor 


analysis. 
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Introducción 


El Perú, a nivel mundial, es uno de los países más afectados por la 
variabilidad y cambio climático (Vega, Lavado, € Felipe, 2018), que 
repercuten en la distribución espacial de los recursos hídricos en la 
cordillera de los Andes de Perú, fuente fundamental para el 
abastecimiento de agua a la población, la agricultura y la producción de 
energía (Wongchuig, Mello, 8 Chou, 2018). En este contexto, se genera 


la necesidad de un manejo sostenible para las cuencas hidrográficas 
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(Bhattacharya, Chatterjee, € Das, 2019), que va desde una 
caracterización morfométrica eficiente (Ghosh 8 Gope, 2021) con nuevas 
herramientas científicas para permitir la adaptación a la variabilidad y el 
cambio climático global con atención en la predicción hasta la solución de 
problemas (Jacobs € Brian, 2020). 


El método de análisis de componentes principales (PCA) permite 
dentro de las diversas variables morfométricas explicar la parte más 
significativa de la varianza de los datos (Helness, Damman, Sivertsen, 8 
Ugarelli, 2019) para identificar relaciones causa-efecto; zonas de clima 
homogéneo o características hidrológicas similares. Este enfoque analítico 
agrupa datos similares mientras identifica relaciones entre variables 
(Gorgoglione, Gioia, € lacobellis, 2019) precisamente para identificar las 
más dominantes en los procesos hidrológicos (Balbín et a/., 2020). 
Subyani, Qari y Matsah (2012) realizaron PCA y análisis de clúster (CA) 
pero solo a 10 cuencas hidrográficas de diferentes tamaños en la región 
occidental de Arabia Saudita con 18 variables morfométricas. Yunus 
Oguchi y Hayakawa (2014) presentaron la cuantificación geomorfológica 
del terreno montañoso, pero limitada a 36 cuencas hidrográficas en la 
península arábiga occidental. Existen otras investigaciones, como la de 
Gajbhiye y Sharma (2017), que solo se basan en la priorización de cuenca 
en función del PCA y el factor compuesto (CF), este último como el cálculo 
del promedio aritmético simple de los rangos de prioridad final, pero 
tienen como desventaja que asumen el mismo grado de importancia a 
todos los parámetros morfométricos (Aher, Adinarayana, 8 Gorantiwar, 
2014; Bharath et al., 2021). El aspecto antes mencionado resulta 
contradictorio, porque cada cuenca hidrográfica tiene sus propias 


características de forma, relieve y red hidrográfica. Investigaciones 
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recientes de Malik et al. (2019), y Setiawan y Nandini (2021) han 
demostrado un optimismo estratégico entre las técnicas PCA y WSA para 
la priorización de cuencas. El método WSA es conocido porque ofrece 
consistencia en el análisis de superficies terrestres que contribuye a la 
priorización de cuencas. Singh y Singh (2018) enfatizan que el WSA es 
un método tradicional de priorización de cuencas en comparación con el 
resto, y brinda resultados dinámicos y eficientes. Nuevos métodos de 
integración son esenciales y permiten una mejor gestión de las cuencas 
hidrográficas en el Perú debido a su diversidad climática y frágiles 
ecosistemas. Por lo tanto, el objetivo de la investigación es priorizar 
cuencas hidrográficas para la conservación de suelos y aguas con base en 
técnicas SIG, PCA y WSA en diferentes regiones del Perú. 


Materiales y métodos 


Se seleccionó como área de estudio a 14 regiones con diferentes 
condiciones hidrográficas, hidrológicas, geológicas y ecológicas que 
integran las 91 cuencas. Dichas cuencas hidrográficas están ubicadas en 
dos de las cuencas hidrográficas más grandes del país: la del Pacífico y la 
del Atlántico. Las cuencas digitalizadas se procesaron en el SIG del 
software ArcMap 10.5, y con base en las fórmulas de la Tabla 1 se pueden 
calcular los parámetros morfométricos de las cuencas. La caracterización 
de las cuencas hidrográficas se realizó en términos cualitativos y 
cuantitativos, que se agrupan en parámetros de superficie, forma, drenaje 


y relieve. 
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Tabla 1. Fórmulas de los parámetros morfométricos de las cuencas 


Parámetros Ecuación Unidad Referencias 
Amplitud de la cuenca W = AL¿* km Horton (1932) 
Pendiente de la cuenca Se = ed nan % Zavoianu (1985) 

Coeficiente de compacidad 
K. = 0.282P40* - Zavoianu (1985) 
de Gravelius 
Factor de forma Ff = AL, - Horton (1932) 
Índice de forma Sy =12A71 % Horton (1945) 
Relación de elongación Re = 1.1284A*"L* - Schumm (1956) 

, Miller, Ritter y Kochel 

Relación de circularidad R¿ = 4TrAP7? - 
(1990) 
Relación de alivio R, = HL;* - Schumm (1956) 
Pendiente media de la 
Sy =4HI5* % Zavoianu (1985) 
corriente 
Densidad de drenaje Di =L ¡45 km/km2 Horton (1945) 
Longitud media del flujo 
Lp =05D7" km Horton (1945) 
terrestre 

Constante de 

E=1D7P km/km2 Schumm (1956) 
mantenimiento del canal 
Coeficiente de rugosidad C, = DaS, - Horton (1945) 


Análisis de componentes principales 
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Es un proceso difícil y complicado clasificar áreas potenciales con más 
beneficios de los recursos hídricos en 91 cuencas. El objetivo del análisis 
de componentes principales (PCA) es explicar la estructura de varianza- 
covarianza en múltiples conjuntos de datos utilizando algunas 
combinaciones lineales de las variables originales, según Kottegoda y 
Rosso (2008). La técnica PCA está diseñada para transformar p variables 
X correlacionadas, que son conocidas, en un número igual de índices Z no 
correlacionados. PCA es una técnica para resumir un conjunto complejo 
de datos, distinguir variables cuantitativas dependientes e 
independientes, así como identificar variables estrechamente 
relacionadas. La notación matricial en su forma general se representa 
mediante la Ecuación (1), en la que Z y X son matrices de nxp y A es una 


matriz de coeficientes de pxp: 
£j = Xaj para j¡=1,2,...p (1) 


Donde Z; es un vector nx1 (columna) y aj es un vector px1 


(columna) de coeficientes., j de 1, 2,... p. 


Considerando que cada columna es una variable observada en rigor 
de n respuestas a la misma pregunta, se pueden estudiar la varianza y 
covarianza entre distintas variables, representadas por una matriz 
cuadrada C de orden pxp simétrica (Yang et al., 2020). La matriz de 
varianza-covarianza X se representa mediante las ecuaciones (2) y (3). 


En la práctica se estima mediante la matriz de covarianza muestral C: 
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Crd. a Cip 

E O C2p 

a : , : (2) 
Cp1 Cp2 Cop 

sl 
Cij¡ = 7 2k=1XijX (3) 


Donde C se convierte en la matriz de correlación de la muestra; X” 
es la transposición de las variables X estandarizadas, y c;j; es el valor de 


cada elemento de la covarianza. 


En el PCA se maximiza la varianza del componente principal (PC), 
que se relaciona principalmente con los valores propios de la matriz de 
covarianza o correlación muestral (Malik et a/., 2019; Kottegoda € Rosso 
2008). La correlación r entre el Xi; de la PC se obtiene por el cociente de 
la covarianza de la matriz X y la varianza de Z; de las PC usando la 


Ecuación (4): 


Cov(x;;z;) 


r(x; z;) = vary (4) 


Donde r es el coeficiente de correlación; Cov (x;, z¡) es la covarianza 
de las variables X en la ¡-ésima fila y la j-ésima fila de los PC; y Var (2;) 
es la varianza de las j-ésimas PC. Un valor de r superior a 0.90 indica una 
fuerte correlación; r entre 0.75 y 0.90, una buena correlación, y r entre 
0.60 y 0.75 señala una correlación moderada (Sharma, Gajbhiye, € 
Tignath, 2015). 


2023, Instituto Mexicano de Tecnología del Agua. | Tecnología y ciencias del agua, ISSN 2007-2422, 


Open Access bajo la licencia CC BY-NC-SA 4.0 14(6), 40-75. DOI: 10.24850/j-tyca-14-06-02 
(https: //creativecommons.org/licenses/by-nc-sa/4.0/) 


a 6) Check for updates 
OPEN ACCESS 


Tecnología y 


CienciaszAgua 


Factor análisis 


El análisis factorial (FA) es un procedimiento similar al análisis de 
componentes principales para identificar variables físicamente 
significativas y reducir parámetros morfométricos. El análisis mediante 
esta técnica produce resultados fácilmente interpretables. Las técnicas de 
FA están bien descritas en la literatura (Finch, 2020). La correlación entre 
las variables está relacionada con las cargas factoriales. Considerando 
correspondencia entre la matriz de puntuación factorial y los Zs en los 
PCs, la ortogonal en la rotación tiene impacto en las PCs y los factores 


correspondientes son equivalentes. 


Por lo tanto, es posible una comparación con los componentes no 
rotados, mientras que la rotación debería ser tal que las interpretaciones 
físicas significativas sean posibles a partir de los componentes resultantes 
obtenidos (Jolliffe, 2002). Por lo tanto, los componentes rotados se dan 


usando la Ecuación (5): 
Z, =XS“1RPD 02 (5) 


Donde Z, son los componentes rotados; D-%-3, una matriz diagonal 
en la que los elementos distintos de O son las raíces cuadradas de los 
valores propios de la matriz de correlación S expresada por la Ecuación 


(2), y RP es la rotación ortogonal. 
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Análisis de conglomerados 


Se aplicó la técnica de análisis de conglomerados jerárquica (CA) para 
clasificar los grupos homogéneos denominados conglomerados. El método 
aplicado para la aglomeración de datos fue en modo R y el criterio de 
Ward (Ward, 1963), uno de los más utilizados en hidrología y 
meteorología (Cupak, Watega, 8 Bogustaw, 2017). El objetivo del método 
fue lograr similitud entre cuencas a partir de criterios morfométricos 


parámetros con enfoque en recurso agua y suelo. 


Análisis de suma ponderada 


WSA se aplicó a las variables morfométricas más significativas que se 
identifican con PCA, WSMP o el peso de los parámetros morfométricos 
que se utilizó en el análisis de correlación cruzada (Siddiqui, Said, Y 
Shakeel, 2020); mientras que el factor compuesto (CF) se evaluó para la 
clasificación de prioridad final, los datos relacionados y categoría, que se 
expresan según Singh y Singh (2018), Malik et a/. (2019), y Setiawan y 
Nandini (2021) como: 


_ Sumof correlation coef ficient 
Wsup UN ; (6) 
Total of correlation 
CF = PRsmpWsmp (7) 
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Donde CF es el factor compuesto; PRsmp es el rango preliminar del 
parámetro morfométrico más significativo del PCA, y Wsmp es el peso del 


parámetro morfométrico significativo. 


Con base en el valor CF, a la cuenca se le asignó una clasificación 
de prioridad. El valor CF más bajo pertenece al rango de prioridad uno, el 
segundo valor más bajo al rango de prioridad dos y así sucesivamente 
para todas las cuencas. El estudio utilizó el paquete estadístico SPSS 26.0 


y STATISTICA 12 como herramienta computacional para el PCA. 


Resultados y discusión 


Los resultados del análisis estadístico en la Tabla 2 de los parámetros 
morfométricos de forma, como el coeficiente de compacidad, factor de 
forma, relación de elongación y factor de circularidad indican la existencia 
de cuencas alargadas y circulares para sus valores medios entre 0.31 y 
1.42, prevaleciendo las cuencas redondas ovales según a los criterios de 
Horton (1932) y Schumm (1956). En el aspecto geológico, el factor de 
circularidad osciló entre 0.24 y 0.73; los resultados altos son indicativos 
de estructuras rocosas que son impermeables y no homogéneas, y que 


controlan el drenaje. 


Tabla 2. Análisis estadístico de parámetros morfométricos. 


Parámetros Mín Máx Media Sd Cv 


Zmáx 124,45 5 500.00 3 698.89 1724.57 46.62 
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Parámetros Mín Máx Media Sd Cv 
Zmín 4.89 4 013.00 1 516.20 1 183.44 78.05 
Hm 45.56 5 429.55 2 182.67 1.395,39 63.93 
A 45.90 1 496.85 459.69 372.18 80.96 
P 33.02 258.12 101.92 52.75 51.75 
Lo 20.40 268.38 85.02 12.60 35.98 
L; 49.53 1 467.02 449.55 364.45 81.07 
Se 6.90 50.99 23.34 13.70 58.69 
W 2.22 21.89 11.30 5.57 49.26 
Kc 1,17 2.06 1.42 0.18 12,49 
Fr 0.10 0.46 0.31 0.09 28.32 
Sw 2.17 9.90 3.68 1.74 47.47 
Re 0.36 0.77 0.62 0.01 15.81 
Re 0.24 0-73 0.51 0.11 21.09 
Ra 0.82 189.18 66.99 44.38 66.25 
Sr 2.36 3/.31 111.63 8.90 76.53 
Da 0.98 3.05 1.88 0.39 20.81 
Lf 0.16 0.31 0.28 0.07 23108 
E 0.33 1.02 0.56 0.13 23.46 
Cr 1.49 143.78 23.91 25.89 108.28 


El relieve accidentado se evidencia con pendientes promedio de 


23.34 %, clasificadas como cuencas montañosas propensas a procesos 
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ambientales dominantes como el deterioro de la fertilidad del suelo, la 
erosión y el transporte de sedimentos en los cursos de agua (Kumar, Das, 
8 Das, 2020). Según Amare, Kassie y Sulla (2020), variaciones 
significativas en el relieve influyen en la distribución de las lluvias y otros 
factores climáticos. La densidad de drenaje en las cuencas hidrográficas 
varía de 0.98 km/km? a 3.05 km/km? para un promedio de 1.88 km/km?. 
Según la categorización de Horton (1945), las cuencas tienen una 
tendencia de drenaje medio y en consecuencia una permeabilidad e 
infiltración de agua media, así como tendencias moderadas a avenidas 
que inducen inundaciones; resultados que están en correspondencia con 
los obtenidos por Mahala (2019), pero en las cuencas hidrográficas del río 


Kosi del norte de la India para el ambiente tropical plano montañoso. 


En general, se trata de cuencas hidrográficas con mayor 
probabilidad de aprovechamiento espacial y temporal de los recursos 
hídricos; por lo tanto, existe una alta necesidad de obras hidráulicas para 
mitigar los problemas de inundaciones y frecuentes derrumbes en las 
cuencas hidrográficas del territorio peruano. La Tabla 3 y Tabla 4 
muestran los resultados del rango preliminar (PRsmp) con base en 


parámetros morfométricos por región. 


2023, Instituto Mexicano de Tecnología del Agua. | Tecnología y ciencias del agua, ISSN 2007-2422, 


Open Access bajo la licencia CC BY-NC-SA 4.0 14(6), 40-75. DOI: 10.24850/j-tyca-14-06-02 
(https: //creativecommons.org/licenses/by-nc-sa/4.0/) 


a y) Check for updates 
OPEN ACCESS 
Tecnología y 


CienciaszAgua 


Tabla 3. Rango preliminar de prioridad basado en las cuencas. 


9 a o 15 S o E AS 
L <= (19) 5 4) o + = 
mn S 1818 56 (8 (5/3 |(£|5 [0 |(8 [8 | $3 
E N [9 [(E|E (8 |[8 |5 |S|2|£ £ 8 | £ 
o ec I£ (a | |O |£ 3 5 | 3 “3 ad [e |ÑS 
o < << S ss ¡5 q 7 
a z] 0) 
E 
Zmáx 2 13 7 9 11 8 14 6 4 10 | 12 1 3 5 


L, 2 8 |[13| 7 12 5 4 14 | 3 1 9 4 6 6 


W 1 5 12| 8 14 3 11 | 13| 4 2 9 6 7 10 


Sc 11 9 8 2 7 12| 14 | 13 | 3 1 5 4 6 10 


Sw 14/11 4 7 1 13| 3 6 8 12 | 10 2 9 5 


Re 11/10 2 7 12 6 9 3 13 5 1 14 8 4 


Lf 6 3 4 | 12 1 2 |14 | 11 9 8 7 5 13 | 10 
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Tabla 4. Clasificación preliminar de prioridad basada en regiones. 


Regiones CP PRsmpP 
Amazonas 5.53 2 
Ancash 7.79 8 
Apurímac 8.37 10 
Cajamarca 7.89 9 
Cusco 8.74 12 
Huancavelica 6.89 5 
Huánuco 9,21 14 
Junín 8.68 11 
La Libertad 6.37 4 
Lima 5.42 1 
Loreto 8.79 13 
Pasco 5.89 3 
San Martín 6.58 6 
Ucayali 7.78 7 


El PCA se inició obteniendo la matriz del coeficiente de correlación 
de Pearson lineal simple (r); existe una fuerte correlación (r > 0.9) entre 
los parámetros morfométricos de 1) A, P; A, Lc; A, Lr y A, W; 2) P, Lc; P, 
L, y P, W; 3) Lc, Lr y Le, W; 4) Lr y W; 5) Kc y Re; 6) Fr, Sw y Fr, Re; 7) Sw; 
Re; 8) Da, Lr y Da, C, y 9) Lr y Cc. Mientras que existe una buena correlación 
(0.75 < r < 0.9) entre 1) Hm, Ra; y 2) Sr, Cr. Además, existen coeficientes 
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de correlación de Pearson moderadamente aceptables (0.60 < r < 0.75) 
en 1) Zmáx, Zmin; Zmáx, Ra; 2) P, Kc; P, Re; 3) W, Fr; W, Sw y W, Re. Los 
parámetros de S¿ y S- no mostraron relación con el resto a las demás 
variables que explican su efecto, con excepción de C,. Resultados 
similares lograron Yunus et al. (2014) en cuencas hidrográficas de la 
península arábiga occidental, y Mokarram y Sathyamoorthy (2015) en 
cuencas hidrográficas de Irán, fundamentalmente en los parámetros de 
A, Le, Re, Fr, Re y Da con fuerte correlación. Mientras que con cierta 
aproximación por los de Malik et al. (2019) está fuertemente 
correlacionado con Ls, Fr, Re, Kc, Re y difiere de Dy en las cuencas de la 
India. Con base en las correlaciones significativas no se observan 
correlaciones de drenaje con las variables superficiales, las cuales 
concuerdan con las alcanzadas por Sharma et al. (2015). Otros autores, 
como Gede, Anwar y Lasminto (2017), demostraron que A y L+, influyen 
en el caudal máximo y el hidrograma unitario en su momento pico con un 
coeficiente de correlación de 0.98, aspecto que podría estar relacionado 


en la presente investigación. 


En la aplicación de la FA a los 20 parámetros morfométricos en la 
primera iteración se definieron seis PC con autovalores mayores a uno y 
que representan alrededor del 88.26 % de la varianza total. Sin embargo, 
los parámetros morfométricos con cargas débiles (Sc y Sr) fueron 
extraídos para mejorar la FA. En estas condiciones se alcanzó la misma 
carga factorial con 6 PC que explican el 94.58 % de la varianza total. La 
rotación Varimax mostró mejora en las cargas (Tabla 6), lo cual se debe 


a que se consideró la matriz de la PC rotada. 
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En la Tabla 5, la carga factorial de PC1 constituye el 26.56 % de la 
varianza total en la matriz rotada; la segunda PC2, el 17.87 %, y la 
tercera PC3 explica el 16.82 % de la varianza total, teniendo estos tres 
primeros componentes una carga fuerte del 60.95 %; mientras que de 
PC4 a PC6, las cargas factoriales solo explican el 33.63 % de la varianza 
total. 


Tabla 5. Matriz de carga factorial de variables morfométricas reducidas 


(Varimax rotado). 


Pmstes Componente principal 

Pci PC2 PC3 PC4 PC5 PC6 

Zmáx 0.15 0.09 -0.10 0.75 -0.03 0.63 
Zmín -0.01 0.04 -0,13 -0.05 -0.07 0.98 
Hm 0.18 0.08 0,01 0.97 0.02 -0.05 
A 0.93 0.26 -0.03 0.01 0.21 0.01 
P 0.83 0.32 -0.02 0.02 0.43 -0.01 
Lo 0.94 0.18 -0.03 0.01 0.21 0.01 
L- 0.93 0.27 -0.04 0.01 0.20 -0.01 
wW 0.83 0.52 0.01 0.02 0.17 0.01 
Kc 0.32 0.01 -0.05 -0.01 0.93 -0.06 
Fr 0.28 0.94 0.06 0.03 0.01 0.02 
Sw -0.26 -0.93 -0.05 -0.06 -0.01 -0.05 
Re 0.29 0.95 0.06 0.03 0.01 0.02 
Re -0.28 0.01 0.04 0.07 -0.92 0.03 
Ra -0.28 -0.02 0.03 0.93 -0.09 -0.04 
Da 0.02 -0.08 -0.97 -0.01 0.08 0.07 
Lr 0.01 0.04 0.99 -0.03 -0.01 -0.05 
C -0.01 0.03 0.98 -0.03 -0.01 -0.05 
Cr -0.50 0.04 LS 0.34 0.28 -0.14 
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La carga factorial de PCi en la Tabla 5 está fuertemente 
correlacionada con las cargas factoriales de A, P, Lc, Lr y W, que podrían 
llamarse factores geométricos que explican el 26.56 % de la varianza total 
de los 18 parámetros analizados. El segundo PC2 tiene una fuerte 
correlación con las variables de Ff, Sw Y Re, por lo que se le denomina 
factor de forma de las cuencas hidrográficas, que explica el 17.87 % de 
la varianza total; mientras que el factor de carga PC3, su fuerte 
correlación es con el Da, Lf y C, denominado factor de drenaje, que explica 
el 16.52 % de la varianza total. La cuarta carga factorial de PC4 tiene una 
fuerte correlación con las variables de Zmáx y Hm, por lo que se denomina 
factor de relieve de las cuencas, que explica el 13.82 % de la varianza 
total. Los primeros cuatro componentes esclarecen el 74.77 % de la 
varianza total, lo que explica los parámetros morfométricos más 
importantes y los factores físicos más significativos. Siendo de gran 
utilidad el PCA para definir las variables más relevantes. Yunus et al. 
(2014) lograron resultados similares a los actuales en 36 cuencas 
hidrográficas, en las que tres PC explicaron el 73 % de la varianza total, 
que reflejaba fuertemente las dimensiones y superficie de la cuenca 
hidrográfica, así como la textura del drenaje. En ese mismo sentido, 
Sharma et al. (2015), tres PC definieron el drenaje, la pendiente y la 
forma; posteriormente, con resultados similares de Malik et a/. (2019), 


con fuerte correlación en forma y parámetros de drenaje. 


La Tabla 6 muestra que existen variables con una fuerte correlación 
en cada PC. Estos parámetros se consideran los más significativos para la 


aplicación del modelo WSA y a su vez la priorización de las cuencas 
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hidrográficas. En la priorización final de las cuencas hidrográficas por 
región, el valor compuesto se calculó con base en el rango preliminar y el 
peso de las seis variables Lc, Re, Lf, Hm, Kc Y Zmín con el análisis de 
correlación cruzada (Tabla 6). El valor compuesto fue determinado 


usando la Ecuación (8). 
CF = 0.26L¿ + 0.23R¿ + 0.101 + 0.16H + 0.18K¿ + 0.08Z pin (8) 


Tabla 6. La correlación cruzada entre los parámetros importantes. 


Parámetros Lc Re Lf H Kc Zmín 
Lc 1.00 0.43 -0.01 0.17 0.48 0.02 

Re 0.43 1.00 0.10 0.16 0.11 0.04 

L; -0.01 0.10 1.00 -0.03 -0.07 ULLZ 

Am 0.17 0.16 -0.03 1.00 0.09 -0.11 

Kc 0.48 0.11 -0.07 0.09 1.00 -0.13 
Zmín 0.02 0.04 -0.17 -0.11 -0.13 1.00 
Sum 2.08 1.84 0.82 1.28 1.47 0.63 
Sum total 8.12 8.12 8.12 8.12 8.12 8.12 
WSA 0.26 0.23 0.10 0.16 0.18 0.08 
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Priorización de cuencas utilizando PCA-WSA 


La clasificación de prioridad final de las cuencas hidrográficas fue de 
acuerdo con el valor de CF más bajo, como se muestra en la Tabla 7 y la 
Tabla 8, donde se establecen cinco clasificaciones de categorías de 
prioridad. La región Amazonas, con categoría de prioridad muy alta de < 
3.10 y área de 11.12 %; la región Lima, con categoría de prioridad alta 
de 3.11 < 5.10, pero con el área mínima afectada de 8.37 %; las regiones 
de La Libertad, Huancavelica, Pasco, Ancash y San Martín, con categoría 
de prioridad media 5.11 < 7.10 y área de 29.71 %; las regiones de 
Apurímac, Junín, Ucayali, Huánuco y Loreto, con categoría de muy baja 


prioridad y el mayor porcentaje de 39.35 % en las menos afectadas. 


Tabla 7. Clasificación prioritaria final de las cuencas hidrográficas en 


función del valor del FC en las regiones. 


Regiones CF Rango de prioridad final 
Amazonas 3.07 1 
Ancash 6.22 6 
Apurímac 9.15 10 
Cajamarca 8.62 9 
Cusco 8.28 8 
Huancavelica 5.58 4 
Huánuco 9.83 13 
Junín 9.33 11 
La Libertad 5.57 3 
Lima 5.07 2 
Loreto 10.45 14 
Pasco 5.59 5 
San Martín 7.09 7 
Ucayali 9.55 12 
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Tabla 8. Categoría de prioridad para las cuencas por regiones. 


, Nivel de Categoría de a Area 
Núm. E EOS Regiones 

prioridad prioridad (%) 

1 < 3.10 Prioridad muy alta Amazonas 11.12 

2 3.11 a5.10 Alta prioridad Lima 8.37 


La Libertad, Huancavelica, 
3 5.11 a 7.10 Prioridad media , 29.71 
Pasco, Ancash, San Martín 


4 7.11a9.10 Baja prioridad Cusco, Cajamarca 11.45 


Apurímac, Junín, Ucayali, 
5 > 9.10 Muy baja prioridad , 39.35 
Huanuco, Loreto 


El mapa de clasificación de prioridades para las cuencas 
hidrográficas por regiones del Perú se muestra en la Figura 1. En la 
investigación, el nivel de categoría más alto representa el grado más alto 
de potencial de erosión y áreas potenciales para la aplicación de la 
conservación de suelos y aguas, así como una necesidad de conservación 
inmediata a la repoblación forestal. En este sentido, las regiones de Lima 
y Amazonas se encuentran en la categoría de alta y muy alta prioridad, 
con un área totalmente vulnerable del 19.49 %. Según Puno y Puno 
(2019), y Setiawan y Nandini (2021), se deben aplicar técnicas 
estructurales de conservación de suelo y agua en las cuencas más 


vulnerables para minimizar la susceptibilidad de las mismas. 
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Figura 1. Mapa de clasificación de prioridades para las cuencas por 


regiones del Perú. 
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La categoría de prioridad media y baja es indicativa de un grado 
moderado de erosión y existencia parcial de cobertura vegetal en siete 
regiones del Perú: La Libertad, Huancavelica, Pasco, Ancash, San Martín, 
Cusco y Cajamarca, representando un área del 41.16 %. En estas cuencas 
se deben implementar programas de reforestación, y un mejor uso 


sostenible de los recursos suelo y agua. 


En la categoría de prioridad muy baja se indica que las cuencas 
hidrográficas tienen características morfométricas sostenibles, como las 
regiones de Apurímac, Junín, Ucayali, Huánuco y Loreto; en esta última 


se debe mantener la protección de la cobertura vegetal del suelo. 


Características hidrológicas homogéneas 


La Figura 2 muestra el dendrograma obtenido del análisis de 
conglomerados en modo R; todos los parámetros de la cuenca se 
presentan con la exclusión de cargas ligeras (S. y Sr). Se observa que 
existen dos conglomerados principales (A y B) y cada uno se subdivide en 
dos conglomerados principales. El conglomerado A se subdivide en el 
conglomerado A(I), con dos parámetros de cuencas (Fr y Re), y el 
conglomerado A (II) incluye seis parámetros (A, Lr, Lc, P, W y Ko); 
mientras que el clúster B se subdivide en el clúster B (III), que incluye 
dos parámetros (Lf y C); y el clúster B (IV), que contiene ocho parámetros 
(Zmáx, Zmin, Hm, Re, Ra, Sw, Da y Cr), el último, con la correlación más baja 
y la distancia más alta. En este sentido, el clúster A presenta mejor 


correlación y menor distancia por debajo de 0.2 con respecto al clúster B. 
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Según Eltahan, Elhamid y Abdelaziz (2021), aunque FA no es un método 
de agrupación, existen similitudes entre el clúster principal A (1 y II) y CP 
(1 y 2), así como B (I y II) y CP (3, 4, 5 y 6). 


5 
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Figura 2. Resultados del análisis de conglomerados en modo R. 
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En la Figura 3 se muestra el resultado del dendrograma obtenido 
por el CA para las 91 cuencas hidrográficas del Perú, numeradas en el eje 
de abscisas. Se observa que existen dos conglomerados principales: C y 
D. El conglomerado C se subdivide en C-I y C-II, con 19 y 13 cuencas, 
respectivamente. El conglomerado D también se subdividió en el 
conglomerado D-III y D-IV en las cuencas 11 y 48, respectivamente. Del 
CA de la Figura 2 y Figura 3, se puede observar que el 35 % de las cuencas 
hidrográficas presenta áreas homogéneas en factor geométrico y de 
forma, de las cuales el 72 % pertenece a la vertiente Atlántico. Sin 
embargo, el 65 % son menos homogéneos en función del factor de 
drenaje y relieve, lo que implica regiones con un 25 %, que tienen 
caudales constantes en los ríos, producto de la alta tendencia a la 
escorrentía en las cuencas ubicadas en las regiones de Loreto, San Martín 


y Cajamarca. 
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Figura 3. Dendrograma obtenido por CA para las 91 cuencas de drenaje 


utilizando Ward Linkage. 


El drenaje y el relieve son la principal causa de mayor disimilitud 
entre las cuencas estudiadas. Esta anomalía se debe a que las cuencas 
hidrográficas ubicadas en la vertiente atlántica y en la llanura amazónica 
tienen una mayor densidad de drenaje que la vertiente Pacífica. Sin 
embargo, la zona norte del Perú debido al fenómeno El Niño Oscilación 
del Sur (ENOS) se ve favorecida por las intensas lluvias, dando lugar a 
mayores crecidas en los ríos, en relación con la zona sur, donde son más 


escasas. En estudios previos, Subyani et al. (2012) y Eltahan et a/. (2021) 
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han obtenido en el CA una similitud de las cuencas hidrográficas en las 


variables geométricas. 


Conclusiones 


La priorización de las cuencas hidrográficas es un factor elemental para 
el desarrollo sostenible y la adecuada gestión de los recursos hídricos. En 
esta investigación se utilizaron técnicas SIG, PCA y WSA para caracterizar 
los parámetros morfométricos de las cuencas hidrográficas, con un 
enfoque en la conservación del suelo y el agua. El estudio reveló que las 
cuencas hidrográficas analizadas en las regiones Lima y Amazonas son 
las más vulnerables a la erosión y presentan áreas potenciales para la 
aplicación de conservación del suelo y agua. Estos resultados contribuyen 
a la toma de decisiones y recomendaciones prioritarias para la 


planificación y gestión integral de las cuencas hidrográficas del Perú. 


La integración de los métodos GIS, PCA y WSA fue exitosa. La WSA 
permitió definir los pesos para los parámetros más relevantes. Asimismo, 
con el PCA se logró determinar la mayor cantidad de parámetros 


significativos, como Lc, Re, Lr, Hm, Kc Y Zmín- 


Las características hidrológicas homogéneas por el dendograma 
explicaron que solo el 35 % de las cuencas hidrográficas son homogéneas 
con base en parámetros geométricos y en forma, con disimilitudes en las 
fuentes de recursos hídricos. El estudio brinda información y herramientas 
útiles para el desarrollo de proyectos de manejo integrado de cuencas 


fluviales a nivel nacional. Asimismo, proporciona una serie de datos 
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importantes para la ejecución de obras hidráulicas para la protección y 


control de ríos en regiones específicas. 


Abreviaturas 


% = Porcentaje. 
AH = Altura media de la corriente. 
A = Área de la cuenca. 


C = Constante de mantenimiento del canal. 


CA = Análisis del clúster. 

CF = Factor compuesto. 

Cov = Covarianza del componente principal. 
Cr = Coeficiente de rugosidad. 

Cv = Coeficiente de variación. 

D-0.5 = Matriz diagonal. 

Dd = Densidad de drenaje. 

e = Equidistancia entre curvas de nivel. 
ENSO = El Niño Oscilación del Sur. 

FA = Análisis factorial. 

Ff = Factor de forma. 

GIS = Sistema de información geográfico. 
Hm = Amplitud altimétrica. 


K¿ = Coeficiente compacidad de Gravelius. 
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km = Kilómetro. 
Km? = Kilómetros cuadrados. 
L = Longitud de la cuenca. 
Lf = Longitud media de la corriente. 
l; = Longitud del contorno. 
L. = Longitud de la corriente. 
Máx = Máximo. 
Mín = Mínimo. 
P = Perímetro de la cuenca. 
PC = Componente principal. 
PCA = Análisis de componentes principales. 
PR = Rango preliminar. 
r = Coeficiente de correlación. 


R2? = Coeficiente de determinación. 


Ra = Relación de alivio. 
Re = Relación de circularidad. 
Re = Relación de elongación. 


RP = Rotación ortogonal. 
S = Varianza. 


S2 = Covarianza. 


Sc = Pendiente de la cuenca. 


Sa = Desviación estándar. 
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S, = Pendiente media de la corriente. 
Sw = Índice de forma. 
W = Amplitud de la cuenca. 
WSA = Análisis de la suma ponderada. 
X = Matriz. 
X' = Variable traspuesta. 
Z = Vector. 
Zmáx = Altitud máxima. 
Zmín = Altitud mínima. 


Z, = Componente rotado. 
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